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Preface 

This  report  describes  the  construction  of  a  miniature  aerodynamic 
skin  friction  drag  balance  and  its  use  in  a  preliminary  investigation 
of  the  drag  reducing  effect  of  boundary  layer  irradiation. 

The  balance  is  of  the  null-position  type  and  is  operated  using 
an  automatic  control  system  to  maintain  the  null  position  of  the  drag 
plate.  High  accuracy  and  stability  are  observed  at  flow  velocities  up 
to  0.8  Mach  Number,  and  a  resolution  of  less  than  0.5%  of  changes  in 
drag  force  can  be  achieved. 

In  preliminary  experiments  irradiation  of  the  boundary  layer  with 

147 

8  particles  from  a  120  mCi  Pm  source  has  given  drag  reductions  of 
5  ±  3%  at  flow  velocities  up  to  0.6  Mach  Number. 

In  future  work  more  accurate  data  on  the  drag  reduction  will  be 
obtained  and  the  physical  basis  of  the  phenomenon  studied. 


1  Introduction 


The  problem  of  reducing  the  drag  of  a  body,  such  as  an  aircraft, 
moving  through  a  fluid  can  be  attacked  either  by  improving  the  stream¬ 
lining  of  the  body,  i.e.  reducing  the  "pressure"  or  "form"  drag  (this 
will  not  be  further  discussed)  or  by  altering  the  fluid/surface  interaction 
in  such  a  way  as  to  reduce  the  local  shear  stress  ("skin  friction")-  at 
the  surface*’^. 

The  region  of  the  fluid  in  which  shear  stresses  become  important 
is  generally  restricted  to  a  thin  layer  on  the  surface  of  the  body,  the 
boundary  layer,  in  which  the  fluid  velocity  rises  from  zero  at  the 
surface  ("no  slip"  condition)  to  the  free  stream  value.  (For  a  detailed 
discussion  of  boundary  layer  theory  see,  for  example.  Refs.  1,  3,  4.) 

It  follows  that  the  properties  of  this  boundary  layer  and  the  adjacent 
surface  are  of  paramount  importance  in  determining  the  drag  forces  acting 
on  the  body  and  are  the  focus  of  any  investigation  into  skin  friction 
drag  reduction. 

In  this  work  an  investigation  has  been  made  of  the  effect  of 
irradiation  of  the  boundary  layer  by  ionising  particles.  The  incentive 
for  this  new  approach  comes  from  work  performed  using  a  torsion  disc 
viscometer  in  an  investigation  of  the  effect  of  ionisation  on  the 
viscosity  of  a  number  of  gases^’^.  (Note  that  in  these  references  the 
activity  level  given  for  the  Po  coated  disc  is  incorrect;  the  lower 
value  of  0.5  to  1.0  mCi  cm-2  is  probably  nearer  to  the  correct  activity.) 
Some  significant  changes  in  the  viscosity  were  observed,  the  most 
prominent  being  an  increase  of  5.2%  for  Ar  at  1  atm  and  a  7.6%  reduction 
for  dry  (i.e.  'v  5  vppm  H^O)  air  at  1  mm  Hg  pressure. 

The  present  research  is  investigating  this  effect  under  conditions 
which  more  nearly  simulate  those  of  flight,  i.e.  gas  static  pressure  close 


to  1  atmosphere  and  flow  velocities  up  to  'v  0.8  Mach  Number  with  both 
moist  and  dry  air  in  addition  to  other  gases. 

To  investigate  methods  of  drag  reduction,  one  must  be  able  to 
measure  skin  friction  drag  accurately.  There  are  several  possible 
approaches  to  this  notoriously  difficult  problem  and  we  briefly  review 
the  most  important.  (For  a  more  detailed  discussion  and  an  account 
of  a  direct  measurement  of  skin  friction  see  Ref.  7.) 

For  flow  satisfying  the  "no  slip"  condition,  the  skin  friction 
shear  stress  is  given  by 


This  expression  is  valid  for  both  laminar  and  (due  to  the  presence  of 
the  sub-boundary  layer)  turbulent  flows.  It  follows  that  a  determination 
of  viscosity  and  of  the  velocity  profile  in  the  boundary  layer  will  yield 
a  value  for  the  skin  friction.  However,  it  is  difficult  (particularly 
if  the  boundary  layer  is  thin)  to  accurately  determine  the  velocity 
profile  and  errors  in  3u/3y  may  be  appreciable.  (It  should  be  noted 
that  most  instruments  commonly  used  for  velocity  measurements,  for 
exanqple,  the  pitot  tube  and  the  hot-wire  anemometer,  are,  due  to  their 
size,  of  limited  use  close  to  the  surface.) 

Heat  transfer  measurements  can  also  (by  virtue  of  the  Reynolds 
analogy  between  heat  and  momentum  transfer)  be  used  to  measure  shear 
stress '*7,  as  with  hot-film  shear  gauges  for  example.  These  have  the 
drawback  of  requiring  extensive  calibration  and  would  not,  in  any  case, 
be  appropriate  to  an  investigation  in  which  the  shear  stresses  on 
different  types  of  surface  were  required. 

This  brings  us  to  the  most  straightforward  method  (in  principle) 


of  skin  friction  measurement,  that  of  a  direct  measurement  of  the  shear 
stress  acting  on  a  small  moveable  element  (drag  plate)  of  the  aerodynamic 
surface;  the  drag  force  being  determined  by  the  application  of  a  measured 
opposing  force  to  the  moveable  element.  The  system  must  be  extremely 
sensitive  (typical  shear  stresses  for  laminar  subsonic  flow  might  range 
from  0.01  -  1.0  gf  cm-2)  and  this  can  make  it  difficult  to  achieve  high 
accuracy. 

An  important  consideration  with  this  method  of  measurement  is  the 
gap  around  the  moveable  element;  if  it  results  in  a  flow  perturbation, 
as  would  be  the  case  when  a  pressure  difference  exists  between  opposite 
sides  of  the  element,  inaccurate  results  will  be  obtained.  The 
alignment  of  the  drag  plate  will  obviously  be  an  extremely  critical 
factor;  the  face  of  the  plate  must  be  accurately  flush  with  the  surface. 

A  misalignment  of  as  little  as  10  pm  could  change  drag  measurements  by 
several  percent.  For  a  detailed  quantitative  discussion  of  these  types 
of  error  sources  see  Ref.  8.  With  regard  to  the  way  in  which  the  force 
is  balanced,  it  is  possible  either  to  measure  the  deflection  of  the  plate 
against  some  resilient  system  (spring  flexures  or  torsion  wires  for 
example)  or  to  maintain  the  plate  at  a  null  position  by  the  application 
of  a  measured  force  (e.g.  magnets  and  coils).  The  null  position  method 
has  the  advantage  that  the  gap  geometry  is  unaltered  and  that  slow  changes 
in  the  elastic  properties  of  supporting  and  resilient  members  are  not 
important. 

Other  factors  which  must  be  considered  are  the  effects  of  temperature, 
external  vibration  and  flow  fluctuations  on  the  calibration  and  response 
of  the  balance  system.  The  best  system  for  this  work  seems  to  be  a  null 
position  drag  balance  mounted  in  the  test  section  of  a  miniature 
windtunnel  designed  to  give  accurately  reproducible  flows,  and  this 
approach  has,  in  fact,  been  used.  For  simplicity,  the  moving  element 


5 


or  drag  plate  forms  part  of  the  wall  of  the  windtunnel  used.  A  more 
detailed  description  of  the  balance  and  tunnel  is  given  in  Section  2. 
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Description  of  Drag  Balance  System 

The  drag  balance  system  consists  of  a  miniature  blow-down 
windtunnel  with  a  sensitive  skin  friction  balance  of  the  null  position 
type  mounted  in  the  test  section.  The  moveable  drag  element  forms  part 
of  the  windtunnel  wall.  The  small  size  of  the  tunnel  makes  it  easy  to 
construct  and  permits  the  use  of  different  working  gases  with  compressed 
gas  cylinders  as  the  supply. 

The  emphasis  has  been  on  developing  a  balance  which  can  reliably 
detect  small  changes  in  the  surface  shear  stress  on  different  types  of 
drag  element  over  a  wide  range  of  flow  velocities,  rather  than  on  an 
absolute  determination  of  the  skin  friction  coefficient.  This  enables 
the  balance  and  windtunnel  to  be  of  a  particularly  simple  form  as, 
provided  the  flow  can  be  accurately  reproduced,  a  detailed  knowledge 

* 

of  boundary  layer  shape,  side  wall  effects,  pressure  gradient  effects 
and  so  on  is  not  essential. 

With  the  present  system,  the  resolution  of  a  change  in  drag  of  as 
little  as  is  thought  to  be  possible. 

By  making  some  approximations,  it  is  possible  to  derive  a  theoretical 
relationship  between  the  measured  drag  and  the  Reynolds  Number  of  the 
flow;  in  the  range  of  validity  of  the  approximations  good  agreement 
is  reached  with  the  drag  data  obtained  in  this  work. 


i.e.  pressure  gradient  effects  as  they  relate  to  the  airflow  pattern. 

A  possible  source  of  error  with  this  design  is  the  pressure  difference 
acting  on  the  ends  of  the  drag  plate.  Given  the  small  pressure  gradient 
along  the  plate,  the  small  pressure  difference  between  the  flow  and 
the  housing,  and  the  gap  size,  relatively  large  in  this  context,  it 
was  felt  that  no  significant  pressure  difference  would  be  maintained 
over  an  appreciable  area  of  the  plates  ends.  In  the  near  future,  the 
magnitude  of  this  effect  will  be  measured  directly. 
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2. 1  Wind tunnel 

The  miniature  windtunnel  used  in  this  work 
is  intended  to  give  two-dimensional  isentropic 

laminar  flow  at  the  test  section  with  low  levels  of  turbulence  and  a 
static  pressure  near  to  atmospheric  for  flow  velocities  of  up  to 
0.8  Mach  Number.  ‘ 

For  simplicity,  the  tunnel  is  of  the  "blow-down"  type  and  runs  with 
either  air  from  a  compressor  which  provides  up  to  68  cu.  ft. /min. 

(2.5  m3/min)  at  approx.  100  psi  (700  kPa) ,  or  from  compressed  gas 
cylinders  containing  pure  gases  or  dry  (zero  grade)  air  at  2000  psi 
(1.4  MPa)  and  having  a  volume  at  s.t.p.  of  6.38  m3. 

The  essential  features  of  the  system  are  shown  in  fig.  1.  The 
pressurised  gas  is  fed  via  appropriate  pressure  regulators  (usually  two 
in  series  to  reduce  the  effect  of  inlet  pressure  changes  on  flow  rate) 
into  the  turbulence  suppressor  tank  (maximum  pressure  10  psi  or 
70  kPa).  This  cylindrical  copper  tank  contains  a  "honeycomb"  array 
of  over  1,400  thin  walled  glass  tubes  of  approximately  0.46  m  in  length 
and  also  two  screens  of  100  mesh  stainless  steel  cloth  at  the  outlet 
end.  This  combination  should  remove  any  large  scale  turbulence  from 
the  airflow.  A  fairing  inside  the  tank  then  guides  the  air  into  the 
rectangular  stagnation  section. 

The  cross-sectional  area  of  the  flow  is  then  reduced  further  to 
that  of  the  balance  section  in  a  two  stage  transition,  the  final  cross- 
section  being  largely  determined  by  the  volume  flow  rate  obtainable  from 
the  gas  source  and  the  desire  for  two  dimensional  flow  over  as  wide  a 
section  of  the  wall  as  possible. 

The  wall  shape  in  each  transition  stage  is  two  circular  arcs  of 
opposite  curvature  (see  fig.  2).  This  shape  is  easy  to  machine  and 
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gives  a  smooch  reduction  in  channel  width. 

The  reduction  from  51  mm  x  126  mm  at  the  stagnation  section  to 
approximately  5.0  mm  x  20.0  mm  for  the  first  balance  and 
approximately  5.0  mm  x  30.6  mm  for  the  second  balance  is  achieved  over 
a  distance  of  about  0.3  m.  (The  increase  in  width  of  the  second 
section  reduces  the  maximum  speed  somewhat,  but  decreases  the  importance 
of  side  wall  boundary  layers  at  lower  speeds.)  This  section  of  the 
tunnel  is  constructed  from  polymethylmethacrylate  (PMMA). 

After  the  balance  section  comes  a  flexible  stainless  steel  duct 
of  0.15  m  diameter  which  exhausts  to  the  outside  of  the  building.  This 
permits  the  use  of  asphyxiating  or  explosive  gases,  and  prevents  any 
hazard  arising  in  the  unlikely  event  of  leakage  from  the  radioactive 
sources  used  in  the  experiments.  (It  should  be  noted  that  all  work, 
even  with  air  and  no  radioactive  sources,  is  conducted  with  the  duct  in 
place,  as  its  presence  alters  the  flow  parameters  sufficiently  for  there 
to  be  apparent  change  in  measured  drag  of  the  same  order  as  the 
resolution  of  the  system.) 


2.2  Prototype  Skin  Friction  Drag  Balance 

A  schematic  diagram  of  the  prototype  null  position  drag  balance  is 
given  in  fig.  3,  and  a  photograph  in  fig.  A. 

The  moveable  element  is  a  steel  plate  of  dimensions  73.0  mm  x 
19.1  mm  x  6.9  mm  fixed  to  a  supporting  steel  plate  (drilled  for 
lightness)  and  suspended  by  flat  steel  springs  attached  to  a  framework 
of  Invar  Steel.  Damping  of  the  plate  is  provided  by  small  glass  dashpots 
containing  an  oil  of  high  viscosity  (normally  used  for  vacuum  work) 


9 


acting  on  an  aluminium  vane  attached  to  the  drag  plate  support. 

The  clearance  between  the  plate  and  the  PMMA  wall  of  the  tunnel 
ia  approximately  0.50  mm  front  and  back  and  0.25  mm  at  the  sides.  The 
plate  is  levelled  by  removing  the  balance  and  top  section  of  the  tunnel 
as  one  piece  and  then  adjusting  three  screws  on  the  supporting  structure, 
alignment  being  checked  with  dial  gauges.  The  balance  housing  is  made 
of  PMMA  and  is  sealed  onto  the  base  by  a  rubber  gasket  which  ensures 
airtightness. 

Position  indication  for  the  plate  is  provided  by  a  Schaevitz  Linear 
Variable  Differential  Transformer  (LVDT)  having  a  range  of  ±  0.010" 

(±  0.25  mm);  readout  is  obtained  from  a  Schaevitz  CAS-025  Signal 
Conditioner  connected  to  a  voltmeter.  The  output  is  roughly  calibrated 
at  15  mV/pm  displacement  by  displacing  the  plate  a  known  distance  using 
feeler  gauges.  The  null  position,  0  V  output  from  the  transformer,  is 
maintained  to  within  'v  5-10  pm  (depending  on  flow  velocity)  by  adjusting 
the  current  in  a  solenoid  coil  (800  turns  of  38  SWG  polyurethane  coated 
copper  wire,  23  fl  resistance)  acting  on  a  J"  (6.35  mm)  AlNiCo-6  permanent 
magnet  attached  to  the  drag  plate.  The  coil  current  is  limited  to 
200  mA  (giving  a  coil  surface  temperature  of  55°C)  corresponding  to 
*  15  gf  drag.  Thermocouples  are  used  to  monitor  coil  and  housing 
temperatures. 

The  average  flow  velocity  at  the  balance  is  determined  from  the 
pressure  at  the  stagnation  section  by  the  usual  expression  for  an 
isentropic  contraction 


P  +  P 


AT 


AT 


Y/(Y  -  I) 


where  P  is  the  gauge  stagnation  pressure  and  the  static  pressure  at  the 
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balance  is  assumed  to  be- atmospheric.  The  stagnation  pressure  is  read 
on  a  precision  manometer  and  can  typically  be  determined  to  an  accuracy 
of  **»  0.2%.  Calibration  of  the  balance  is  performed  using  a  thin 
(0.004"  or  0.1  mm)  flat,  constant  strength  cantilever  spring  with  four 
strain  gauges  cemented  to  it,  giving  a  resistance  bridge  with  four 
active  arms.  The  gauges  are  M-line  MA-06-1 25AC-350  (steel  compensated) 
and  readout  is  from  a  Phillips  PR9307  Carrier  Wave  Measuring  Bridge. 

The  spring  is  itself  calibrated  by  being  clamped  horizontally  and 
different  known  weights  (analytical  balance  weights)  hung  from  it.  The 
calibration  of  the  balance  is  then  achieved  by  moving  the  spring  (attached 
to  the  drag  plate  by  a  nylon  thread)  along  a  microscope  type  slide 
(see  figs  5,  6,  7).  This  approach  was  adopted  in  preference  to  the 
more  obvious  methods  in  the  hope  that  the  absence  of  sliding  or  rolling 
friction  in  the  calibration  system  might  lead  to  a  superior  calibration. 
However,  although  adequate  calibrations  could  be  obtained,  the  non-linearity 
of  the  calibrations,  the  double  calibration  procedure  and  the  drift  of 
the  spring  calibration  turn  out  to  outweigh  the  potential  advantages 
and  this  method  was  not  subsequently  used. 

Other  problems  experienced  with  the  prototype  balance  are  difficulty 
of  control  (to  hold  a  particular  flow  velocity  and  accurately  adjust 
the  current  in  the  coil  at  the  same  time  is  a  difficult  two  man  job); 
the  levelling  procedure  is  not  sufficiently  accurate;  adjusting  the  plate 
to  hang  squarely  in  the  gap  in  the  tunnel  wall  and  without  any  friction 
between  coil  and  magnet  or  LVDT  and  core  is  an  extremely  difficult 
procedure;  some  magnetic  sticking  also  occurred  probably  due  to  slight 
magnetisation  of  the  drag  plate  or  steel  base  of  the  balance.  These 
latter  factors  make  it  hard  to  get  an  adequate  degree  of  reproducibility 
of  results  on  changing  drag  plates.  The  former  problem  has  been  overcome 
by  the  development  of  an  automatic  control  system  (see  Section  2.3). 


II 


Data  taken  with  the  prototype  balance  indicates  that  the  basic 
design  is  satisfactory  and  an  improved  model  has  been  constructed  which 


does  not  have  the  drawbacks  mentioned  above.  This  is  described  in 
Section  2.4. 


2.3  Automatic  Control  System 

As  mentioned  in  Section  2.2, it  is  difficult  to  control  the  proto¬ 
type  balance  by  hand;  the  regulator  requires  constant  fine  adjustment 
to  cope  with  inlet  pressure  changes  and  to  control  the  coil  current 
at  the  same  time  to  an  acceptable  degree  of  accuracy  requires  two  people. 
Accordingly  a  control  circuit  was  designed  which  would 
interface  between  the  LVDT  signal  conditioner  and  the  coil  power  supply 
and  act  in  such  a  way  as  to  maintain  the  null  position  of  the  plate, 
i.e.  OV  output  from  the  LVDT  (see  fig.  8).  The  power  supply  has  a 
high  impedance  input  socket  accepting  voltages  of  *v<  -  I  V  to  +  1  V, 
giving  proportional  output  currents  from  a  power  transistor  of 
^  -  250  mA  to  +  250  mA  with  a  20  0  load.  The  control  circuit  must 
therefore  increase  its  output  voltage  if  the  LVDT  voltage  is  positive, 
decrease  its  output  voltage  if  the  LVDT  voltage  is  negative,  and  hold 
its  output  if  the  LVDT  voltage  is  zero  (or  close  to  zero).*  The 
circuit  must  also  be  stable  with  regard  to  external  excitation  (flow 
fluctuations,  floor  vibrations) . 

It  can  be  seen  that  a  direct  feedback  system  is  not  appropriate 
and  would  make  the  system  unstable.  An  integrating  system  would,  in 

A  downstream  displacement  gives  a  +ve  output  requiring  a  +ve  current 
to  balance  it;  the  null  position  is  set  to  be  close  to  zero  by  sliding 
the  LVDT  in  its  mount;  small  tare  or  zero  currents  through  the  coil 
then  hold  it  at  the  null  position. 
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theory,  work,  but  in  the  absence  of  vibration  might  be  slow  to  approach 
zero;  with  vibration  it  is  possible  (depending  on  the  phase  lag  and 
time  constant  of  the  system)  that  the  system  could  be  unstable. 

A  system  was  designed  which  would  make 
adjustments  in  small  stages,  giving  the  balance  time  to  respond  to 
each  stage  before  proceeding.  A  circuit  building  up/down  a  fine 
staircase  voltage  for  positive/negative  inputs  was  used; 
this  has  the  same  response  speed  even  at  voltages  close  to 
zero.  After  some  initial  design  work, the  circuit  shown  in  figs  9,  10 
was  constructed  and  found  to  perform  satisfactorily. 

The  LVDT  output  is  fed  (via  an  RC  filter  which  removes  the 
2.5  kHz  ripple  on  the  LVDT  output  and  other  high  frequency  signals)  into 
an  amplifier  stage  giving  a  gain  of  up  to  'u  I03.  The  amplifier  output 
stages  are  then  used  to  trigger  Schmitt  Nand  gates  whose  output  gates 
clock  pulses  from  a  timer  chip.  (The  circuit  is  symmetrical  over  this 
stage,  different  halves  of  the  circuit  being  used  to  handle  positive 
and  negative  signals.)  The  pulses  are  then  fed  via  opposed  diode- 
zener  diode  combinations  to  a  capacitor  bank;  the  diodes  prevent  leakage 
of  charge  when  the  gated  clock  outputs  are  at  zero.  The  capacitors  are 
of  the  polycarbonate  type  which  have  a  very  low  leakage  rate  (see  fig. 

10). 

The  circuit  therefore  takes  the  input  signal,  amplifies  it,  and, 
if  large  enough,  triggers  the  passage  of  voltage  pulses  of  the  appropriate 
polarity  to  the  capacitor  bank  where  each  pulse  adds  or  removes  a  step 
from  the  staircase  output  voltage.  For  an  input  between  the  trigger 
levels,  the  output  voltage  remains  almost  constant;  the  decay  rate 
being  negligible  for  output  voltages  in  the  designed  working  range  of 
±  1  V.  Trigger  levels  can  be  adjusted  using  the  10  kfl  and  500  kft 
potentiometers  (zero  offset  and  gain  adjust)  in  the  amplifier  stage  to 
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anything  from  a  few  mV  to  a  few  volts;  for  use  with  the  drag  balance 
they  are  set  at  ±  3.5  mV  (±  II  mV  with  the  prototype).  The  clock  chip 
frequency  and  mark/space  ratio  may  be  adjusted  using  the  2  Mfl 
potentiometers.  For  use  with  the  drag  balance, the  frequency  is  set 
to  70  Hz  with  a  pulse  width  (-ve  at  the  clock  stage)  of  2.7  ms  and 
a  dead  time  of  11.5  ms. 

An  evaluation  of  the  performance  of  the  drag  balance  control  system 
is  given  in  Section  3.2. 


2.4  Drag  Balance 

The  principle  of  the  drag  balance  is  the  same  as  for  the  prototype. 
Photographs  showing  the  windtunnel  system  and  balance  are  given  in 
figs  11-14. 

The  drag  plate  (brass,  74.4  mm  x  19.35  mm  x  6.85  mm)  is  attached 
to  the  "crossbar"  on  an  H  shaped  bakelite  holder  supported  from  a  rigid 
duralumin  box  by  12  wires  (steel  piano  wire  0.010",  0.25  mm,  in  diameter 
and  'v  100  mm  long);  up  and  down  at  the  end  of  each  arm  and  side  to  side 
at  the  centre  of  each  arm.  This  constrains  the  plate  to  motion  parallel 
to  the  flow  direction  only.  The  duralumin  housing  has  ports  for  access 
top  and  bottom,  sealed  with  perspex  windows  and  0  rings,  and  panels  at 
the  front  for  electrical  and  pressure  lead  throughs.  The  supporting 
wires  are  raked  slightly  (^  1°)  and  are  fixed  to  tensioning  screw 
adjusters.  This  permits  extremely  accurate  levelling  and  location  of 
the  drag  plate  within  the  hole  in  the  perspex  wall  of  the  windtunnel; 
the  clearance  is  about  0.25  mm  front  and  back  and  0.10  mm  at  the  sides. 
The  wires  are  adjusted  to  have,  as  far  as  possible,  similar  tensions; 
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this  can  be  achieved  by  tuning  them!  The  tensions  are  such  that  the 
resonant  frequency  of  the  balance  is  'v  30  Hz.  Plate  changing  is 
accompl xshed  by  removing  a  section  of  the  underside  tunnel  wall  and 
dropping  the  drag  plate  through. 

Plate  alignment  is  by  an  electrical  resistance  method,  using  a 
levelling  plate  (see  fig.  15)  consisting  of  four  brass  contacts  set  in 
an  insulating  block.  The  block  has  two  ridges  which  are  held  against 
the  upper  windtunnel  surface  and  the  brass  contacts  (accurately  coplanar 
with  the  ridges)  line  up  with  the  edges  of  the  drag  plate.  A  9  V 
potential  is  applied  between  the  drag  plate  and  the  contacts;  electrical 
contact  is  indicated  by  the  light  emitting  diodes  (LED’s)  which  can  be 
seen  in  the  right  hand  PMMA  panel  on  the  balance.  The  tensions  in 
the  supporting  wires  are  adjusted  until  all  four  LED's  flicker  as  the 
plate  swings.  This  has  proved  to  be  an  extremely  accurate  and  reproducible 
method  of  aligning  the  drag  plate. 

The  LVDT  used  has  a  range  of  ±  0.025"  (0.65  mm)  and  the  bridge 
output  is  calibrated  at  ^  20  mV/ym. 

The  drag  force  is  balanced  by  two  solenoidal  coils  of  1350  turns 
of  38  SWG  polyurethane  coated  copper  wire  on  a  12  mm  o/d,  1  mm  thick, 
aluminium  former  25  mm  long  acting  on  10  mm  diameter  AlNiCo-6  magnets. 

This  allows  much  higher  drags  to  be  balanced,  permitting  investigations 
in  the  turbulent  regime.  The  coils  are  connected  in  parallel  (to  give 
equal  forces)  and  have  resistances  of  56  and  57  0.  The  coils  are 
enclosed  in  aluminium  housings  filled  with  a  heat  sink  compound  and  fixed 
to  the  duralumin  wall  of  the  housing  which  acts  as  a  heat  sink. 
Thermocouples  are  used  to  monitor  coil  temperatures;  compared  with  the 
previous  coil  only  negligible  temperature  rises  are  recorded.  Adjustment 
of  the  clearance  between  the  core  and  casing  of  the  LVDT  and  also  of  the 
magnets  and  coils  is  achieved  by  pivoting  them  about  points  which  are 
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nominally  lined  up;  if  the  alignment  is  not  sufficiently  accurate 
clearance  is  obtained  merely  by  altering  the  angles  at  which  the 
components  are  clamped. 

Damping  is  by  the  same  method  as  before  but  with  a  larger  vane 
area  giving  a  damping  ratio  of  about  0.14. 

The  ACS  maintains  the  plate  position  to  better  than  5  pm  during 
steady  running. 

The  centreline  flow  velocity  at  the  balance  section  is  determined 
from  a  pitot  tube  and  static  orifice  15  mm  downstream  from  the  drag 
plate.  (It  is  assumed  that  there  is  no  transverse  pressure  gradient.) 

The  pitot  tube  has  a  slightly  smoothed  flat  end  and  has  an  o/d  of 
1.25  mm,  and  i/d  of  0.65  mm,  and  a  length  of  'v*  25  mm  before  the  stem. 

The  static  orifice  has  the  same  diameter  as  the  pitot  tube.  The  static 
pressure  at  the  plate  and  the  difference  between  this  and  the  housing 
pressure  are  also  measured;  the  latter  gives  a  sensitive  check  for  leaks 
from  the  housing.  The  flow  static  temperature  and  balance  housing 
temperature  are  measured  using  thermocouples.  Calibration  of  the  balance 
is  performed  by  hanging  weights  from  a  thread  attached  to  the  damping 
vane  and  passing  through  the  balance  housing  wall  and  over  a  pulley 
(see  fig.  16).  The  pulley  is  small  and  light  and  is  mounted  using 
miniature  precision  bearings.  The  vreights  used  are  accurately  weighed 
(0.6970  ±  0.0003  gf  each)  steel  spheres  (small  ball  bearings)  and  are 
handled  using  a  slightly  magnetised  rod  to  prevent  contamination.  The 
use  of  bearings  rather  than  analytical  balance  weights  minimises  the 
disturbance  to  the  drag  balance  during  loading  (up  to  15  gf). 

Extremely  accurate  calibrations  can  be  made  by  this  method,  and  with 
the  new  levelling  and  plate  changing  techniques  a  high  degree  of 
reproducibility  between  data  for  similar  drag  plates  can  be  attained. 


Figure  1  VIEW  OF  WiNDTUNNEL  SYSTEM 


Figure  2 


Side  view  of  wind tunnel 


transition  section 


Schematic  Diagram  of  Prototype 
Skin  Friction  Drag  Balance 
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View  of  Prototype  Skin  Friction 
Drag  Balance 
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Side  View  of  Calibration 
Spring  Assembly 
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Figure  6 


Front  View  of 
Calibration  Spring 
(x  4) 


Schematic  Diagram  of  Calibration 
Procedure  for  Spring 
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Schematic  Diagram  of  Drag  Balance 
Automatic  Control  System 


I. 


SCHEMATIC  DIAGRAMS  OF  CONTROL  CIRCUIT 


CONTROL  CIRCUIT 


Circuit  Diagram  of  Control  Circuit 


Inverting  Operational  Amplifier 

Timer  Chip  in  Astable  Multivibrator  Mode 

Quad  2-Input  Schmitt  Nand  Gates 


Diodes  -  IN4I48  Capacitors  -  polycarbonate 


Figure  1 2 


View  of  Skin  Friction  Drag  Balance 


In  the  left  of  this  picture  can  be  seen  the  second  stage 
of  the  windtunnel  transition  section,  leading  into  the  balance. 

The  two  removable  panels  on  the  front  of  the  balance  are 
for  coil  current  and  thermocouple  lead  throughs  (on  the  left), 
and  for  pressure  and  LVDT  connections  (on  the  right) .  The  top 
of  the  right  hand  panel  holds  the  LED’s  used  in  the  drag  plate 
levelling  process. 

The  knurled  adjusters  (two  on  the  lateral  faces  of  the 
housing,  and  four  top  and  bottom)  are  used  for  tensioning  the 
supporting  wires. 

The  start  of  the  exhaust  duct  and  its  pressure  tap  can  be 
seen  in  the  right  hand  section  of  the  photograph. 


TOP  INSIDE  VIEW  OF  BALANCE 


Figure  13 


Figure  14 
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View  from  below  of  the  inside  of 
the  drag  balance 


The  brass  drag  plate  can  be  seen  through  the  removable  ► 

section  in  the  lower  windtunnel  wall.  Beyond  the  upper  surface 
of  the  tunnel  the  H  shaped  drag  plate  holder  can  be  seen. 
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Photograph  of  Balance  During  Calibration 


The  photograph  illustrates  the  ease  with  which  the  balance 
can  be  calibrated;  the  exhaust  duct  is  removed,  the  pulley  fixed 
in  place,  and  the  plug  in  the  wall  unscrewed.  This  has  the 
cotton  thread  attached  to  the  plate  holder  running  through  it, 
and  the  balance  is  then  ready  for  calibration.  Access  to  the 
interior  of  the  balance  housing  is  not  needed. 
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.3  Experimental  Work 

This  section  describes  the  experimental  work  performed,  consisting 
of  balance  calibrations,  obtaining  drag  curves  for  a  range  of  flow 
speeds,  and  then  taking  drag  data  at  several  sample  flow  velocities 
for  a  range  of  different  drag  plates. 

The  balance  was  recalibrated  and  a  new  drag  curve  obtained  after 
each  plate  change.  Reproducibility  between  similar  plates  is  'v  {%. 
During  the  course  of  the  experiments,  a  check  of  other  experimental 
parameters,  such  as  temperatures  of  housing,  coils  and  flow,  was 
maintained,  but  the  small  range  of  variation  of  these  parameters  seemed 
to  have  no  significant  effect  on  either  the  calibration  or  the  drag 
data. 

In  the  type  of  set  up  used,  it  is  difficult  to  predict  precisely 
the  drag  forces  to  be  expected.  At  high  flow  velocities,  however,  i.e. 
when  the  boundary  layer  is  thin  compared  with  the  channel  thickness,  it 

is  possible  to  approximate  the  situation  by  a  flat  plate  in  a  uniform 

* 

flow  (see  fig.  17  overleaf).  Here  X  is  the  effective  entry  length  of 

the  channel;  one  might  expect  x=0  to  correspond  to  a  point  somewhere 

in  the  transition  section  of  the  windtunnel.  u  will  be  taken  to  be 

00 

equal  to  the  centreline  velocity,  uq,  of  the  tunnel  near  the  drag  plate. 
The  local  shear  stress  on  one  side  of  the  plate  is  given  by* 

T  *=  0.332  /w  u  3/2/^  . 

o 


* 


Assuming  that  the  pressure  gradient  along  the  channel  is  not  too 
great. 


SCHEMATIC  DIAGRAM  SHOWING  FLAT  PLATE  APPROXIMATION 
APPLIED  TO  WINDTUNNEL  AND  DRAG  BALANCE 
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The  drag  on  the  plate  is  then 


D  =  /  idA  =  W  /  Tdx  =  0.664  W  4ip  u  '^(/x  +  1  -  )  . 

plate  x=X  ° 


Thus  the  skin  friction  drag  coefficient  is  given  by 


For  X=0,  this  reduces  to  the  usual  Blasius  equation  for  skin  friction 
drag  on  a  plate  in  uniform  laminar  flow. 

Good  agreement  with  this  formula  has  been  reached  for  6/a  <  0.4; 

X  may  be  determined  from  the  experimental  data  and  as  expected  corresponds 
to  an  entry  length  beginning  near  the  centre  of  the  two  stage  transition. 

At  lower  flow  velocities,  boundary  layers  (including  those  on  the 
side  walls)  are  not  small  compared  with  the  channel  dimensions  and  it 
is  difficult  to  predict  the  relationship  between  and  Re. 


3. 1  Calibration  and  Initial  Tests  of  Prototype  Balance 

Prior  to  calibration  the  strain  gauge  spring  was  loaded  cyclically 
(to  10  gf  in  0.5  gf  steps)  several  times  to  work  out  any  strains  in  the 
gauge  cement.  Once  an  acceptably  low  degree  of  hysteresis  had  been 
achieved  rough  calibrations  were  obtained  by  loading  in  one  direction 
only.  Checks  were  made  on  the  amount  of  drift  under  load;  this  was 
found  to  be  acceptable.  To  estimate  the  effective  tare  weight  of  the 
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pan  and  horizontal  spring,  best-fit  power  curves  were  determined  for  the 
first  six  points  of  each  calibration  curve  for  a  range  of  assumed  tare 
weights.  The  curves  of  r2  vs  weight  gave  a  "pan  load"  of  0.30  ±  0.01  gf. 

A  final  set  of  seven  calibrations  was  then  obtained  (see  fig.  18 
for  a  typical  calibration  curve)  and  best-fit  power  curves  of  the  form 

b 

R/mV  -  as  (1/gf)  8 


calculated. 

Results  were  as  follows:- 


Table  1  Spring  Calibration 


Calibration  Code 

ag 

b 

s 

D1 

0.527 

0.960 

D2 

0.534 

0.971 

D3 

0.530 

0.971 

D4 

0.529 

0.955 

D5 

0.540 

0.966 

D6 

0.530 

0.975 

D7 

0.540 

0.963 

For  each  calibration  r2  >  0.9995. 


There  was  no  evidence  for  an  ag  on  bg  dependence,  a  temperature 
dependence  (within  the  temperature  limits  encountered),  or  a  dependence 
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t 


on  bridge  warm  up  time,  so  final  values  for  a  and  b  are  taken  as 

s  s 

direct  averages  of  the  above  values,  yielding: - 


a  *=  0.533  ±  0.002 

s 

b  =  0.966  ±  0.003 

s 


(standard  errors) 


The  calibration  of  the  balance  is  then  performed  by  establishing 

an  accurately  known  current  in  the  coil  and  bringing  the  LVDT  reading 

to  zero  using  the  calibration  spring,  recording  the  bridge  reading  when 

this  had  been  accomplished.  See  fig.  19  for  a  typical  calibration  curve. 

In  all  ten  accurate  calibrations  were  performed  and  best  fit  power  curves 

bc 

obtained:-  R/mV  =  ac(I/mA) 

Results  were  as  follows:- 


Table  2  Prototype  Balance  Calibration 


Calibration  Code 

a 

b 

c 

c 

El 

0.0382 

0.947 

E2 

0.0504 

0.908 

E3 

0.0399 

0.944 

E4 

0.0425 

0.932 

E5 

0.0507 

0.913 

FI 

0.0371 

0.960 

F2 

0.0406 

0.943 

F3 

0.0426 

0.937 

F4 

0.0522 

0.901 

F5 

0.0408 

0.940 
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Again  r2  >  0.9995  except  for  F3  (0.998). 

These  results  give 

a  =  0.044  with  S.E.  0.002  b  =  0.93  with  S.E.  0.02 

c  c 

The  final  calibration  equation  has  the  form 

1/gf  =  a(I/mA)B 

where 

a 

Error  calculations  for  a  and  6  are  complicated  by  the  dependence  of 
the  second  calibration  on  the  first  and  by  a  correlation  between  a 

c 

and  b  (a  =  0. 31  -  0. 28  b  ) .  Results  for  a  and  6  with  estimated 
c  c  c 

errors  are 

a  »  0.076  ±  0.09  8  ■=  0.96  ±  0.02 

The  calibration  curve  is  shown  in  figure  20.  Some  initial  runs 
with  the  balance  using  zero  grade  air  were  then  made;  the  drag  curves 
obtained  are  shown  in  figures  21  and  22.  At  high  flow  speeds,  good 
agreement  is  reached  with  the  Blasius  type  curve;  X  is  found  to  be 
approximately  0.23  m. 

The  behaviour  of  the  vs  R^  curve  at  R^  106  suggests  the  flow 
is  becoming  turbulent  at  this  point,  but  due  to  the  limitations  of  the 
coil  data  at  higher  flow  could  not  be  obtained. 


l/be 


b  /b 
c  s 


3.2  Performance  Tests  of  Automatic  Control  System 

Charts  recording  plate  position  and  balance  control  circuit 
output  have  been  obtained  for  a  range  of  flow  conditions,  showing 
both  steady  state  behaviour  and  response  of  the  automatic  control 
system  (A.C.S.)  to  changing  flow  conditions,  simulated  by  setting  up 
a  constant  flow  before  switching  the  circuit  on  (see  figs  23,  24,  25). 
The  oscillations  noticeable  in  the  traces  for  high  flow  velocities  are 
a  consequence  of  the  vibration  of  the  steel  exhaust  duct.  Considering 
its  simplicity,  the  control  circuit  exhibits  surprisingly  good  response 
and  stability  characteristics,  and  has  been  found  to  be  entirely 
adequate  over  the  whole  range  of  calibration  and  flow  conditions 
encountered  in  this  work. 


3.3  Drag  Data  Obtained  with  Prototype  Balance  and  ACS 

The  data  shown  in  fig.  26  was  obtained  using  the  ACS  and  "wet" 
air  from  the  compressor.  Data  was  taken  mainly  over  the  range  in 
which  agreement  with  the  Blasius  curve  was  obtained  and  also  around 
Rj  *  I06  to  see  if  a  transition  was  again  observed.  In  fact  no  signs 
of  deviation  from  the  C^  «  curve  were  observed  up  to  the  maximum 

flow  velocity  attained,  ^  0.8  Mach  Number.  It  is  possible  that  the 
previous  data  was  inaccurate  (difficult  to  obtain  good  data  without 
ACS)  but  it  is  more  likely  that  a  transition  to  turbulent  flow  had 
been  taking  place,  resulting  from  either  slight  misalignment  of  the 
drag  plate,  or  misalignment  of  the  test  section  with  the  first 


transition  section. 
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The  good  agreement  between  the  experimental  and  "theoretical" 
curves  and  the  degree  of  reproducibility  of  results  obtained  for  a 
given  plate  configuration  (of  the  order  of  i%)  indicate  that  the 
concept  of  the  drag  balance  is  satisfactory,  and  a  modified  balance 
with  better  reproducibility  between  plates  should  prove  to  be  adequate 
for  drag  reduction  investigations. 


3.4  Calibration  and  Initial  Tests  of  Drag  Balance 

Before  calibrating  the  balance,  it  was  decided  to  determine  the 

resonant  frequency  of  the  balance  and  the  effectiveness  of  the  dashpot 

damping.  To  do  this  a  7.5  mV  rms  signal  from  a  Levell  Oscillator  was 

fed  into  the  input  socket  of  the  coil  power  supply.  The  corresponding 

LVDT  voltage  amplitude  was  measured  on  an  oscilloscope  and  a  plot  of 

amplitude  ratio  for  a  range  of  frequencies  obtained  (see  fig.  27) . 

If  we  call  k/w  the  damping  factor  where  k  and  u>  can  be  defined  from 
o  o 

the  equation  of  motion:- 

X(a)C 

x  +  2kx  +  to  2x  =  le  (the  forcing  function  can  have  unit 

amplitude  w.l.o.g.) 

The  amplitude  ratio  is  then 


A 


1 


and  this  is  a  maximum  for 


to  =  a)  given  by  -gH  =  0 

P  *0) 

P 


w 

=  1  - 


2k2 


to 


The  peak  amplitude  is  then 


A 

P 


_ 1 _ 

(2k/to  )  Jl  -  k2/w  ?- 
o  v  o 


It  follows  that  a  measurement  of  A  and  io  will  give  w  /<o  and  k/to  : 

p  p  p  o  o 


0)  2 

_E_ 

“o2 


F 


n 

p 


Jsl  =  i 
0)  2  2 
o 


to  2^ 

1 - P7 

<0  * 


This  gives 


to  =  29  Hz 

o 


k/to  =  0.14 
o 


(The  theoretical  curve  for  these  values  is  shown  as  the  dashed  curve 
in  fig.  27.) 

The  damping  is  some  way  from  optimum  but  should  be  adequate.  (The 
resonant  frequency  will,  of  course,  change  if  the  tensions  in  the  wires 
are  altered,  e.g.  on  changing  drag  plates,  but  as  the  tensions  are  kept 
roughly  the  same  the  above  values  may  be  taken  as  typical.) 


A  rough  set  of  drag  data  was  then  taken,  including  measurements 
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of  the  static  pressures  at  the  plate  and  pitot  tube  (see  fig.  28). 
It  can  be  seen  that  in  calculating  the  Mach  Number  of  the  flow  from 


P 

PS 

P  +  P 
at  s 

it  is  a  fair  approximation  (errors  a  few  %  only)  to  neglect  Pg  compared 
with  P  (though  it  cannot,  of  course,  be  neglected  compared  with  the 

a  C 

pitot  pressure) . 

During  this  test, the  maximum  leak  pressure  recorded  was  'v  2  mm  Hg 

at  P  g  =  200  mm  Hg.  Removing  the  4  BA  stopper  in  the  hole  in  the  housing 

wall  (for  calibration  purposes)  a  leak  pressure  of  'v  10  mm  Hg  is  obtained. 

There  is  thus  a  sensitive  indication  of  leaks;  a  maximum  leak  pressure 

of  'v-  2  mm  Hg  at  P  =  200  mm  Hg  would  seem  to  be  tolerable, 
ps 

The  plate  was  then  accurately  aligned,  wire  tensions  checked,  and 
a  first  accurate  set  of  data  taken  for  a  range  of  flow  speeds.  During 
this  extended  procedure,  some  drift  of  the  aero  current  is  noticeable. 

If  1^  is  the  current  of  the  n^  reading,  N  the  number  of  readings,  and 
I  and  1^  the  zero  currents  at  the  start  and  finish  respectively  of  the 
test,  the  current  can  be  corrected  by 

!c  -  In  "  IsO  "  n/N)  -  Ifn/N  , 


1  + 


Y  ~  1 


Y  -  1 


-  1 


assuming  the  drift  is  roughly  linear  with  time  as  readings  are  taken. 

The  curve  of  I  vs  P  obtained  is  shown  in  figure  29.  (For  better 
c  ps 

accuracy  in  drag  comparisons,  a  series  of  readings  are  taken  at  a 
particular  flow  velocity,  checking  the  zero  before  and  after  each 
reading.  This  is  done  for  a  sample  of  six  or  seven  flow  velocities.) 
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The  data  shows  that  the  balance  is  working  satisfactorily  and 
so  a  series  of  accurate  calibrations  were  performed  using  bearings  of 
weight  0.6970  ±  0.0003  gf  (the  error  is  the  uncertainty  in  the  average 
weight,  rather  than  the  variation  between  bearings).  A  typical 
calibration  curve  is  shown  in  fig.  30. 

The  results  were  as  follows 


Table  3  First  Balance  Calibration 


Calibration  Code 

Slope  of  Curve 
(mA/bearing) 

r2 

Coil  Constant 
(gf  mA'1) 

C0905802 

2.829 

0.9996 

0.2464 

C0905803 

2.814 

0.9987 

0.2477 

C0905804 

2.838 

0.9998 

0.2456 

C0905805 

2.813 

0.9999 

0.2478 

C0905806 

2.836 

0.9998 

0.2458 

Cl  205801 

2.827 

0.9999 

0.2465 

Mean  Coil  Constant 

=  0.2466  gf/mA 

Standard  Deviation 

=  0.0009  gf/mA 

Standard  Error  =  0.0004  gf/mA 


The  accuracy  achieved  by  this  quick  and  simple  method  is  far  greater 


than  for  the  strain  gauged  spring.  (Some  slight  hysteresis  is  noticeable 
on  sections  of  the  calibration  curves  due  to  the  pulley  not  being 
absolutely  free,  but  the  effect  is  usually  negligible.  Occasionally 


26 


the  pulley  becomes  a  bit  "sticky",  but  cleaning  the  bearings  usually 
rectifies  this.) 

Having  calibrated  the  balance, a  curve  of  drag  coefficient  vs 
Reynolds  Number  could  be  obtained;  this  is  shown  in  fig.  31.  Again, 
agreement  is  achieved  with  the  modified  Blasius  curve  at  high  flow 
velocities.  (The  curve  is  higher  than  that  obtained  with  the  prototype; 
this  is  to  be  expected  as  the  channel  is  wider  and  the  inlet  section 
has  been  modified.)  A  value  for  X  of  ^  0.18  m  is  determined  from  the 
data,  x  =  0  corresponding  to  a  point  just  downstream  of  the  midpoint 
of  the  transition  section. 

Accurate  drag  measurements  were  then  made  at  pitot  static  pressures  of 
5,  10,  25,  50,  100,  150,  200  mm  Hg  corresponding  to  approximate  Mach 
Numbers  of  0.06,  0.14,  0.22,  0.30,  0.42,  0.50,  0.58. 

It  should  be  noted  that  accuracy  in  reading  the  coil  current  is,  at  low 
flow  speeds,  limited  mainly  by  external  vibration,  giving  a  resolution 
in  coil  current  of  ±  0.05  raA,  or,  exceptionally  (on  Sundays  for  example!) 

1  0.025  mA.  A  summary  of  the  data  obtained  is  given  below  (coil  currents 
in  mA) . 


Table  4 


First  Set  of  Drag  Measurements 


Pps/mm  Hg  5 
Calibra¬ 

10 

25 

50 

100 

150 

200 

tion  Code 

D 1205801  0.9 

2.1 

9.9 

19.6 

36.3 

49.7 

61.2 

D1 205802  0.9 

2.1 

9.8 

19.7 

36.5 

50.1 

61.2 

D 1205803  0.9 

2.1 

9.8 

19.8 

36.7 

50.2 

61.9 

D 1205804  1.0 

2.0 

9.8 

19.8 

36.7 

50.4 

61.6 

D 1205805  0.9 

2.2 

9.9 

19.8 

36.6 

50.4 

61 .9 

0.9+0.05  2.110.05  9.810.05  19. 7510.05  36.5510. 10  50. 1510.10  61. 5510.15 
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With  a  coil  constant  of  0. 2466±0.0004  gf/mA,this  gives  the 
following  drag  values:- 


Table  5  First  Set  of  Drag  Values 


P  /mm  Hg  D/gf 

ps 


5 

0.22 

+ 

0.01 

10 

0.52 

+ 

0.01 

25 

2.42 

± 

0.01 

50 

4.87 

± 

0.02 

100 

9.01 

± 

0.03 

150 

12.37 

± 

0.03 

200 

15.18 

± 

0.04 

Except  at  the  lowest  flows  this  should  be  sufficiently  accurate  to 
detect  drag  changes  of  a  fraction  of  a  percent.  As  a  quick  check  on 
reproducibility, the  drag  plate  was  changed  for  a  similar  plain  brass 
plate  and  the  same  procedure  followed.  The  corresponding  I  vs  P^s 
curve  is  shown  in  fig.  32. 


I 
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Four  calibrations  were  then  performed  as  follows 


Table  6  Second  Balance  Calibration 


Calibration  Code 

Slope  of  Curve 
(mA/bearing) 

r2 

Coil  Constant 
(gf  mA-1 ) 

Cl  805801 

2.826 

0.9999 

0.2467 

Cl  805802 

2.820 

0.9999 

0.2472 

C1905801 

2.833 

0.9999 

0.2461 

Cl  905802 

2.839 

0.9999 

0.2455 

Mean  Coil  Constant  =  0.2464  gf/roA 

Standard  Deviation  =  0.0009  gf/mA 

Standard  Error  =  0.0004  gf/mA 

This  is  consistent  with  an  unchanged  calibration. 

A  rough  set  of  drag  measurements  were  then  taken  yielding  the  data 
below: - 

Table  7  Second  Set  of  Drag  Measurements  (D1905801) 


P  /mm  Hg  5 

ps 

10 

25 

50 

100 

150 

200 

I/raA  0.9±0.05 

2.05±0.05 

9.8+0.  1 

19.7+0. 1 

N/A 

48.4+0.5 

G1.8±0.5 
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giving  drag  values  of:- 


Table  8  Second  Set  of  Drag  Values 


/mm  Hg 

D/gf 

5 

0.22  ±  0.01 

10 

0.51  ±  0.01 

25 

2.41  ±  0.02 

50 

4.85  ±  0.03 

100 

N/A 

150 

11.9  ±  0.1 

200 

15.2  ±  0.1 

Except  for  P  ■  150  mm  Hg  this  is  in  good  agreement  with  the  first 
ps 

set  of  values.  This  suggests,  as  expected,  that  the  major  problem 
in  working  with  the  radioactive  sources  will  be  attaining  a  sufficiently 
accurate  drag  plate  surface. 

As  a  preliminary  experiment,  it  was  decided  to  investigate  the 
effect  of  the  sources  when  mounted  on  the  opposite  wall  to  the  drag 
plate.  Although  the  irradiation  of  the  boundary  layer  on  the  plate  will 
be  less,  5  mm  of  air  should  not  cut  the  radiation  level  too  greatly  in 
the  case  of  the  a  and  $  sources.  The  6+,  y  source  will  give  essentially 
the  same  level  of  y  radiation  at  the  plate,  but  the  positrons  will  be 
stopped  within  a  very  short  distance  of  the  source. 


The  surface  roughness  of  the  active  plate  is,  of  course,  less 
critical  in  this  configuration,  so  although  the  activity  at  the  drag 
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plate  may  be  lowered  the  ultimate  resolution  of  the  balance  is  achieved. 

The  construction  of  the  radioactive  drag  plates  is  described  in 
the  next  section,  and  the  results  of  this  initial  investigation  in  the 
subsequent  section. 


3.5  Construction  of  Radioactive  Drag  Plates 


The  ideal  approach  to  the  construction  of  active  drag  plates  would 
involve  taking  a  standard  plate  and  coating  it  with  a  thin  film  of  a 
radioactive  isotope.  Unfortunately,  at  present,  this  is  a  prohibitively 
expensive  procedure.  Accordingly,  standard  "static  eliminator"  radioactive 
foils,  obtained  from  the  Radiochemical  Centre  at  Amershan,  have  been 
incorporated  into  the  drag  plates.  The  specifications  of  the  radioactive 
foils  used  are  given  in  Table  9. 

The  polonium  and  promethium  isotopes  are  sealed  onto  a  substrate  by 
a  thin  (<v  3  pm)  protective  coating  of  silver.  The  cobalt  is  infused 
directly  into  a  copper  strip.  For  use  in  the  drag  balance,  the  foils 
are  incorporated  into  a  drag  plate  of  standard  dimensions  in  such  a  way 

as  to  give  an  accurately  flat,  clean  surface.  The  plate  is  milled  out 

* 

over  an  area  slightly  greater  than  that  of  the  foil  ,  and  the  foil  is 
then  cemented  into  the  plate  using  beeswax  (supplied  from  a  dental 
laboratory).  This  is  a  surprisingly  strong  and  versatile  adhesive;  it 
has  the  advantages  that  when  molten  it  flows  readily  and  is  transparent 
(allowing  a  check  on  the  position  of  the  foil),  when  set  it  is  extremely 

Note  that  the  ^'®Po  foil  was  trimmed  slightly  before  use. 


i 
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hard  and  rigid  (although  somewhat  brittle),  but  may  easily  be  softened 
and  removed  if  the  foil  position  is  unsatisfactory. 

Initially  the  foil,  typically  a.  0.5  mm  thick  and  slightly  curved, 
is  fixed  to  an  accurately  flat  brass  backing  strip,  ^  3  urn  thick,  using 
double-sided  sticky  tape.  The  foil  and  milled  plate  are  then  placed  on 
a  flat  surface,  consisting  of  greaseproof  paper  fixed  to  a  flat  glass 
block  with  double  sided  sticky  tape  and  then  wiped  with  3  in  t  graphited 
oil.  The  oil  soaked  paper  prevents  wax  getting  onto  the  surface  of  the 
foil,  and  permits  easy  release  of  the  plate  from  the  surface  once  the  wax 
has  set. 

The  outer  section  of  the  plate  is  clamped  to  the  block  by  spring 

clips  and  the  foil  is  then  aligned  and  clamped  by  screwing  five  screws 

through  the  back  of  the  milled  plate  and  onto  the  backing  strip.  The 

assembly  is  then  carefully  warmed  with  a  blowtorch  and  molten  beeswax 

run  into  the  plate  through  holes  cut  in  the  back.  (A  diagram  showing 

the  plate  during  the  foil  mounting  procedure  is  given  in  fig.  33.  A 

photograph  showing  the  working  area  used  and  a  plate  during  construction 

can  be  seen  in  fig.  34.  It  should  be  noted  that  all  manipulation  of  the 

foils  is  performed  wearing  gloves  and  glasses,  and  that,  where  possible, 

the  foils  are  handled  with  tongs  behind  a  lead  screen  having  a  lead  glass 

window  in  it.  A  film  monitoring  badge  is  always  worn;  no  dosage  greater 

0.02  rem/month  has  been  recorded.  Particular  care  is  taken  not  to  damage 

9)0 

the  protective  surface  of  the  foils,  especially  the  alpha  foil,  “  Po 
being  extremely  hazardous  if  ingested.) 

The  melting  point  of  the  beeswax  is  %  70°C  and  the  maximum  foil 
temperature  reached  during  mounting  is  probably  ^  I00°C,  which  should 
be  perfectly  safe.  Once  the  wax  has  set  hard  (after  about  3/4  hour) 
the  clamping  screws  are  removed  and  the  drag  plate  released  from  the 
surface  and  wiped  clean. 
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Two  silver  foils  and  the  alpha  and  beta  foils  have  so  far  been 
mounted.  Except  for  the  alpha  foil  the  standard  of  surface  achieved 
has  been  good.  In  the  case  of  the  alpha  foil, the  ends  of  the  foil  have 
been  capped,  making  it  hard  to  get  a  flat  surface,  and  on  releasing 
from  the  surface  one  cap  partly  unfolded.  It  has  not  been  possible  to 
refix  this  satisfactorily,  so  this  plate  has  a  slight  protrusion  at  one 
end.  Figure  35  shows  one  of  the  silver  drag  plates  in  the  holder  used 
for  the  experiment  described  in  the  next  section. 


3.6  Effect  of  Boundary  Layer  Irradiation 

In  this  section, an  account  of  a  preliminary  investigation  into  the 
effect  of  boundary  layer  irradiation  is  given.  It  seems  obvious  that  the 
drag  balance  will  be  less  sensitive  to  surface  imperfections  on  the  wall 
opposite  to  the  drag  plate,  than  to  those  on  the  drag  plate  itself.  It 
follows  that  one  of  the  radioactive  plates  could  be  rapidly  mounted,  to 
a  reasonably  high  accuracy,  in  the  lower  wall  of  the  windtunnel.  The 
possible  loss  in  accuracy  is  compensated  for  by  the  reduced  exposure  to 
the  source.  There  will  also,  of  course,  be  a  drop  in  the  radiation  level 
at  the  drag  plate,  but  for  a  5  mm  thick  airflow  this  should  not,  except 
perhaps  for  the  alpha  source,  be  too  great. 

Accordingly,  as  a  preliminary  experiment,  the  active  drag  plates  have 
been  mounted  flush  with  the  lower  windtunnel  surface,  opposite  to  a  plain 
brass  drag  plate,  and  drag  measurements  made  at  several  flow  speeds. 

Each  plate  is  mounted  in  turn,  and  three  complete  sets  of  drag  measurements 


taken  each  time.  The  results  are  given  in  Table  10. 

Each  experiment  involves  measuring  the  zero  current,  drag  current, 
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Table  10  Drag  Measurements  with  Active  Plates 


Drag  values  in  mA  coil  current.  Calibration  is  approximately  0.246  gf/mA. 


Expt. 

Code 

Drag  Plate 
Type 

PPS=5 

10 

25 

50 

100 

150 

200  mm  Hg 

RD01 

0.6 

1.55 

9.7 

19.2 

33.4 

45.3 

56.0 

RD02 

► 

Agtfl 

0.6 

1.6 

9.8 

19.3 

33.6 

45.8 

53.5 

RD03 

> 

0.6 

1.6 

9.7 

19.3 

33.6 

46.0 

57.0 

RD04 

*> 

0.45 

1.75 

9.5 

19.3 

33.6 

45.1 

55.3 

RD05 

► 

Agtf2 

0.5 

1.75 

9.6 

19.5 

34.0 

45.2 

55.5 

RD06 

-- 

0.5 

1.75 

9.65 

19.5 

34.1 

45.4 

55.4 

RD07 

I 

0.3 

1.4 

8.9 

18.3 

31.7 

42.8 

53.0 

RD08 

1 

► 

I47d 

Pm 

0.35 

1.4 

8.9 

18.3 

31.8 

42.8 

52.5 

RD09 

J 

0.4 

1.4 

8.9 

18.3 

31.7 

42.8 

52.9 

RDIO 

] 

0.5 

1.4 

9.4 

19.0 

32.9 

44.6 

55. 1 

RDI  1 

1 

1 

► 

Agttl 

0.45 

1.45 

9.3 

19.0 

33.0 

44.5 

55.0 

RD 1 2 

J 

0.5 

1.6 

9.4 

19.0 

33.1 

44.6 

54.9 
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zero  current,  drag  current  and  then  zero  again,  for  each  flow  velocity 
and  then  obtaining  the  average  drag  current.  The  measurements  are  taken 
in  quick  succession,  but  changes  in  flow  velocity  are  made  sufficiently 
slowly  for  the  plate  displacement  to  be  minimal,  imposing  no  undue  stresses 


on  the 

suspension 

system. 

The 

mean  drag  values 

for  each 

type  of  plate 

have  been  calculated  and 

are  given  below: 

l  ” " 

i 

Table 

1 1 

Mean  Drag 

Values 

Plate 

Type 

p  =5 
‘PS 

10 

25 

50 

100 

150 

200  mm  Hg 

Ag*l 

0.6 

1.6 

9.7 

19.3 

33.5 

45.7 

55.5 

Ag#2 

0.5 

1.75 

9.6 

19.4 

33.9 

45.2 

55.4 

147., 

Pm 

0.35 

1.4 

8.9 

18.3 

31.7 

42.8 

52.9 

Ag*l 

0.5 

1.5 

9.4 

19.0 

33.0 

44.6 

55.0 

The  data  taken  with  the  silver  plates  gives  an  idea  of  the  spread 
of  values  to  be  expected  due  to  variations  in  the  accuracy  of  the  plate 
surface  and  of  the  positioning  in  the  windtunnel  wall.  If  we  assume 

147 

similar  uncertainties  apply  to  the  results  for  the  Pm  foil  then  values 

147 

for  the  drag  reduction  with  the  Pm  plate  can  be  calculated  together 
with  rough  error  estimates.  These  values  are  summarised  below,  and  will 
be  discussed  in  Section  5.  Note  that  the  values  at  Ppg=200  mm  Hg  seem 
surprisingly  accurate;  this  is  because  the  air  compressor  is  delivering 
almost  its  full  capacity  and  so  the  inlet  pressure  at  the  flow  regulators 
changes  only  slowly. 
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Table  12  Summary  of  Results 


Pps/mm  Hg 

Approx.  Mach 

Drag  Values 

Drag  Values 

%  Decrease  in  Drag 

Number 

with  Ag  Plate 

with  Plate 

with  147Pm  Plate 

5 

0.06 

0.55 

±  0.05 

0.35 

±  0.05 

(40 

+ 

20) 

10 

0.14 

1.6 

±  0.1 

1.4 

±  0.1 

(15 

± 

15) 

25 

0.22 

9.6 

±  0.2 

8.9 

±  0.2 

7 

4 

50 

0.30 

19.2 

±  0.2 

18.3 

±  0.2 

5 

± 

2 

100 

0.42 

33.5 

±  0.5 

31.7 

±  0.5 

5 

± 

3 

150 

0.50 

45.2 

±  0.6 

42.8 

±  0.6 

5 

± 

3 

200 

0.58 

55.3 

±  0.3 

52.9 

±  0.3 

4 

± 

1 

210 

A  few  rough  results  taken  with  the  Po  foil,  which  has  a  poor 
surface  with  a  slight  protrusion  at  one  end,  gave  increases  in  apparent 
drag  of  •v  2%,  of  the  order  of  the  uncertainty  of  the  experiment,  and 
little  significance  is  attached  to  this  result. 


Typical  Calibration  Curve  for 
the  Strain  Gauge  Spring 


Typical  Calibration  Curve  of  Prototype  Drag  Balance 
using  Strain  Gauged  Spring 


100  150  200 

COIL  CURRENT  FOR  RESTORING  FORCE/mA 


Final  Calibration  of  Prototype  Drag  Balance 


100  150  200 

COIL  CURRENT  FOR  RESTORING  FORCE/mA 

DRAG  BALANCE  CALIBRATION 


Shear  Stress  vs  Flow  Velocity 


Preliminary  Test  using  Zero  Grade  Air 


ERODYNAMIC  SHEAR  STRESS, /(g force  cm 


Figure  22 

Drag  Coefficient  vs  Reynolds  Number 


as  fig.  21 


Control  Circuit  Output  for  a  Range  of  Trigger  Conditions 


Inputs  causing  either  positive,  negative  or  no  triggering 
were  fed  to  the  control  circuit  to  simulate  the  application  of 
simple  forces  to  the  drag  plate. 

The  output  voltage  of  the  control  circuit  was  recorded  as  a 
function  of  time. 

The  circuit  was  connected  to  the  coil  power  supply  during 
this  experiment. 

The  staircase  structure  of  the  output  is,  of  course,  too 
fine  to  be  resolved  in  this  way.  (Use  of  an  X-Y  plotter  would 
seem  to  be  an  obvious  way  of  recording  this,  but  the  excessively 
low  input  impedance  of  the  plotter  causes  misleading  results.) 


CONTROL  :  CIRCUIT 


Figure  24 

Response  Curves  for  the 
Drag  Balance  Control  System 


RESPONSE  CURVES  FOR  DRAG  BALANCE  CONTROL  SYSTEM 


Performance  Curves  for  the  Automatic  Control  System 
Under  "Steady"  Flow  Conditions 
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Drag  Coefficient  vs  Reynolds  Number 


(Reynolds  Number  is  based  on 
the  plate  length) 
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Figure  27 


Response  Curve  of  Drag  Balance 


tude  Ratio  vs  Frequency  of  Forcing  Function 


I  Graph  Data  Ref.  5943  L09  4  Cycles  *  3  Cycles 


Coil  Current  vs  Pitot  Static  Pressure 


First  Set  of  Test  Data 


Typical  Calibration  Curve  for  the  Drag  Balance 


Coil  Current  vs  Load  (number  of  bearings) 


Drag  Coefficient  vs  Reynolds  Number 


First  Set  of  Test  Data 


wwiwaw’-  b 


Coil  Current  vs  Pitot  Static  Pressure 


Second  Set  of  Test  Data 
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Diagram  Showing  the  Foil  Mounting  Procedure 


Note  that  the  diagram  is  not  to  scale,  and  that  the  spring 


clips  have  been  omitted  for  the  sake  of  clarity 


Photograph  Illustrating  the  Precautions  Taken  When 
Handling  the  Radioactive  Sources 


The  shield  which  can  be  seen  in  the  righthand  side  of  the 
picture  is  constructed  of  2"  thick  interlocking  lead  bricks. 
Behind  the  screen  is  a  drag  plate  clamped  as  for  mounting. 

In  the  left  of  the  picture  can  be  seer,  the  radiation  monitors 
used  to  check  the  degree  of  exposure  to  the  sources,  and  to 
check  for  radioactive  contamination. 


Active  Drag  Plate 


This  photograph  shows  a  drag  plate  made  up  with  a  silver 
foil.  The  plate  is  fixed  into  a  perspex  holder  by  double  sided 
sticky  tape,  ready  for  use  in  the  experiment  described  in 
Section  3.6. 

The  surface  scratches  noticeable  in  the  photograph  are 
in  fact  very  fine,  but  have  been  enhanced  in  the  picture  by  the 
illumination  used  during  the  exposure. 
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4  Summary  of  Results 


The  drag  balance  can  be  readily  calibrated  to  an  accuracy  approaching 
0.1%  and  yields  accurate  drag  measurements  over  the  large  part  of  the 
subsonic  flow  regime  ('v  0. 1-0.8  Mach  Number).  Above  <\»  0.35  Mach  Number 
good  agreement  can  be  obtained  with  the  Blasius  Equation  for  the  drag 
coefficient. 

During  operation  the  position  of  the  drag  plate  is  maintained  to 
within  'v  2  pm  of  the  null  position  and  the  average  position  to  considerably 
better  than  this. 

Drag  data  taken  for  different  plain  drag  plates  shows  reproducibility 
to  a  few  percent. 

In  a  preliminary  experiment, skin  friction  drag  reductions  were 

observed  on  irradiating  the  flow  over  a  plain  brass  plate  with  0.2  MeV 

147 

8  particles  from  a  120  mCi  Pm  source.  The  results  are  summarised 
below.  Note  that  little  importance  should  be  attached  to  the  two  low 
flow  measurements. 


Flow  Mach 
Number 


%  Drag  Decrease 


(40  ±  20) 
(15  ±  15) 
7^4 
5  4  2 


5  ±  3 


5  ±  3 


4  ±  I 


5  Discussion  and  Conclusions 


The  drag  balance  shows  good  accuracy  and  stability  and  the  prospects 
of  achieving  a  sensitivity  of  ±  0.5Z  in  drag  measurements  are  excellent 
except  at  the  lowest  flow  speeds  where  building  vibrations  are  a  problem. 

The  results  obtained  with  @  particle  irradiation  show  that  a 
decrease  in  drag  force  is  taking  place  at  atmospheric  pressure. 

Variations  in  measurements  due  to  differences  in  the  foil  mounting  and 
plate  positioning  will  be  considered  in  more  detail  in  future  research; 
however,  it  is  felt  that  the  silver  and  promethium  plates  are  sufficiently 
alike  for  the  drag  change  observed  to  be  taken  as  a  real  effect. 

Neglecting  the  very  uncertain  data  at  low  velocities  the  results  show 
a  reduction  in  drag  on  irradiation  of  5  ±  3Z.  It  is  tempting  to  interpret 
the  data  as  showing  that  the  change  in  drag  is  smaller  at  higher  flow 
speeds  than  lower.  This  would  agree  with  our  ideas  on  the  mechanism  of 
the  effect  (discussed  below). 

The  changes  in  drag  (or  effective  viscosity)  are  unlikely  to  be 
caused  by  the  introduction  of  Coulomb  forces  due  to  ionisation  in  the 
air;  the  magnitude  of  this  effect  should  be  orders  of  magnitude  smaller 
than  that  observed. 

In  both  this  work  and  that  of  Clark,  Kestin  and  Shah,  a  conducting 
drag  element  was  used;  the  interaction  between  the  ions  and  image 
currents  could  act  so  as  to  reduce  the  viscous  shear,  but  the  magnitude 
of  this  effect  is  not  known. 

A  conjecture  is  that  the  clustering  of  neutral  species  around  the 
ions  might  be  important;  this  obviously  causes  mobility  changes,  and 
it  is  possible  that  the  clusters  might  also  act  in  such  a  way  as  to 
inhibit  the  formation  of  turbulence  at  the  molecular  level.  However, 
at  present,  there  is  insufficient  knowledge  about  the  density,  structure 


and  stability  of  the  clusters  formed  for  any  estimate  to  be  made  of  the 
magnitude  of  this  effect.  Studies  of  the  basic  mechanisms  involved  in 
this  phenomenon  will  form  a  major  part  of  future  research. 
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Appendix  A  Problems  Encountered 

The  work  was  delayed  earlier  this  year  by  the  late  arrival  and 
inoperable  condition  of  the  air  compressor.  Both  the  Phillips  PR9307 
and  the  Schaevitz  CAS-025  bridges  broke  down  during  the  year,  the  repair 
of  the  latter  being  hindered  by  the  unavailability  of  the  correct  circuit 
diagram.  Except  for  these  occasions  the  equipment  has  performed 
adequately. 

It  is  probable  that  a  different  type  of  alpha  foil  will  have  to  be 
used  in  future,  as  it  is  difficult  to  get  a  flat  drag  plate  using  the 
polonium  foil. 
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